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I. INTRODUCTION 

Frictional motion plays a central role in diverse systems 
and phenomena that span vast ranges of scales, from the 
nanometer con t acts in herent in micro- and nanomachines 
(lUrbakh et all. I2004T) and biological molecular motors 



Bormuth et all . l2009| ) to the geophy sical scales charac- 
teristic of earthquakes (Scholz, 1998). Due to its enor- 



mous practical and technological importance, the prob- 
lem has stimulated progress over the centuries. Historical 
figures from Leonardo da Vinci onwards have brought 
friction into the field of physics, with the formulation 
of time-honored phenomcnological frictional laws, which 
have been referred to as the Coulomb- Amontons laws. 
These statements can be briefly summarized as follows: 
(i) frictional force is independent of the apparent area of 
contact; (ii) frictional force is proportional to the nor- 
mal load; (iii) kinetic friction (the force to keep relative 
motion at constant speed) does not depend on the slid- 
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ing velocity and is smaller than static friction (the force 
needed to initiate motion between two contacting bodies 
at rest). Also in the light of a mass of empirical data, 
serious attempts were made in the first half of the 20th 
century toward a microscop ic understanding of these laws 
( Bowden and Taborl . Il950[ ). Whereas the basic physics 
underlying sliding friction - non equilibrium statistical 
mechanics of solids, sheared fluids, and moving surfaces 
- is in principle quite exciting, the field as a whole has 
(even if with notable exceptions) failed to attract ade- 
quate interest by the physicist until the last few decades. 
A lack of microscopic data, and a corresponding lack of 
theory, have perhaps contributed to project an unattrac- 
tive image of sliding friction. 

Three quiet revolutions, of broad nature and unre- 
lated to friction, are radically changing this state of af- 
fairs. First, progress in the general area of complex- 
ity provided new tools to tackle non-equilibrium disor- 
dered systems with interacting degrees of freedom. Sec- 
ond, and crucial, the developments in nanotcchnology 
extended the study of friction and permitted its anal- 
ysis on well-characterized materials and surfaces at the 
nano- and microscale. Notably the invention of scanning 
tip inst ruments of the Ato mic Force Microscope (AFM) 
family (jBinnig et al.l . ll986f) has opened nanofriction 
brand new a venue; the use of t he Surface Force Appa- 
ratus (SFA) (|lsraelachvilil Il992h has led to the system- 
atic studies of confined mesoscopic systems under shear; 
while instru ments such as the Quartz Crystal Mic robal- 
ance (QCM) (jKriml . Il996i: iKrim and Widoml . Il988l) mea- 
sure the inertial sliding friction of adsorbate submono- 
layers. Thanks to these methods, a mass of fresh data 
and information on well defined systems, surfaces, mate- 
rials, ph ysical conditions has accumula ted in the last two 
decades ( Carpick and Salmeronl[l997t ). The resulting in- 
sights into the atomic size contacts themselves in terms 
of chemical interactions and of the elementary processes 
that are involved in the excitation and dissipation of en- 
ergy are changing our perspective. Third, computer sim- 
ulations have had a strong boost, also allowed by the fan- 
tastic growth of computer power. The numerical study 
of frictional models on one hand, and direct atomistic 
molecular dynamics simulations on the other hand, are 
jointly advancing our theoretical understanding. Invalu- 
able initial reviews of the progress brought about by these 
revolutions in our physical underst anding o f sliding fric- 
tion can be found in the books by iPersso u ( 2000a|) and 
iMatd (120081 ). 

Despite the practical and fundamental importance of 
friction and the growing efforts in the field, many key 
aspects of dynamics of friction are not yet well under- 
stood. Even for the most studied nanoscale systems, 
such as AFM sliding on graphite or NaCl surfaces, a mi- 
croscopic mechanism of friction is still lacking, and ex- 
perimental observations (for instance, velocity and tem- 
perature dependencies of friction) have been rational- 
ized within simplified models including empirical param 
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eters (jBarel et a/.l . l201ll;|Jansen et all |2010| ; iRiedo et all 



21 )().'{ 1. Fundamental theory is still difficult in all fields 
of sliding friction, including nanofriction, since the slid- 
ing motion generally involves sudden nonlinear stick-slip 
events, that cannot be treated within traditional theo- 
retical approaches such as linear-response theory and hy- 
drodynamics. Experiments in tribology have long suf- 
fered from the inability to directly observe what takes 
place at a sliding interface. Although AFM, SFA and 
QCM techniques have identified many friction phenom- 
ena on the nanoscale, many interpretative pitfalls still 
result from indirect or ex situ characterization of contact 
surfaces. In this Colloquium we will briefly cover some 
aspects, progress and problems in the current modeling 
and simulation of sliding friction, from nano to mesoscalc. 
In nanoscale friction we consider systems that are small 
enough to be treated at the atomistic scale, such as in 
the AFM experiments. For larger systems we need a 
mesoscopic approach that lies in between the atomistic 
details and the macroscopic behavior. In the spirit of a 
Colloquium we intend to draw examples from our own 
experience to illustrate some of the concepts and points 
of interest, growth, and doubt in selected forefront areas. 

One of the main difficulties in understanding and pre- 
dicting frictional response is the intrinsic complexity of 
highly non-equilibrium processes going on in any tri- 
bological contact, which include detachment and reat- 
tachment of multiple microscopic junctions (bonds) be- 
tween t he surfaces in relative motion while still in 
contact (iBormuth et al. l. l2009HGerde and Marderl l200lt 
lUrbakh et all 12004 ) Therefore friction is intimately re- 
lated to instabilities that occur on a local microscopic 
scale, inducing an occasional fast motion of the corre- 
sponding degrees of freedom even if the slider's ccntcr- 
of-mass velocity is extremely small. Understanding the 
physical nature of these instabilities is crucial for the elu- 
cidation of the mechanism of friction, as we will empha- 
size below. 

Sliding friction has been addressed following different 
types of theoretical approaches: "minimalistic" models 
(MM) atomistic molecular dynamics (MD) simulations, 
mesoscopic multicontact models and phenomenological 
rate-state (RS) models. 

MMs are discussed in Sec. [TTJ They provide an inter- 
mediate level description between atomic scale physics 
and macroscopic phenomenological approaches like RS 
models, focusing on a small number of relevant degrees 
of freedom which describe the sliding motion, and ex- 
hibit instabilities during stick-slip. Applications of MMs 
provided explanations for phenomena of high complexi ty 
[see, e.g., ([Miiser et aUl2003tlVanossi and Braunl 120071 )] ■ 
On the whole, MMs are playing a major role in ratio- 
nalizing the wealth of nano and mesoscale friction data 
produced over the last decades. 

Atomistic MD simulations, discussed in Sec. IIII1 
have a wide range of applicability in nanoscale fric- 
tion, and have reached a high l evel of rigor and accu- 
racy (jRobbins and Mtiserj . l200lh . But, as discussed in 
Sec. IIII.B1 they are mostly limited to time and length 
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scales which are too short to emulate many tribological 
phenomena. An important issue, therefore, is how to re- 
duce the large-scale, many-parameter MD equations to 
simpler mesoscale descriptions based on fewer degrees of 
freedom. 

Multicontact models, discussed in Sec. HVl provide 
such a simplified description in terms of dynamical for- 
mation and rupture of elastically coupled contacts. At 
the largest macroscopic scale, phenomenological RS mod- 
els simplify the description even further, introducing one 
or two dynamical equations with coefficients chosen to 
fit experimental quantities and then used to describe 
a wide range of observed frictional behavior, such as 
the transition b etween stick-slip (regular or chaotic) and 
smooth sliding ( Carlson and Batistal Il996l ). a nd varia- 
tions of friction for a sudden change of velocity ( Dietrich! 
119791: lRuinal . ll983l) . RS models are often the best avail- 
able approaches to describe macrosc opic friction in the 
ordin ary world, from the micro size (iBaumberger et al\ . 
Il999h to larger and larger scales (|Scholzlll998l ). However, 
most of the "state variables" in RS models cannot be eas- 
ily related to physical system properties, a fact that limits 
the insight and predictive power of these models. This 
Colloquium is limited to nano and mesoscale frictional 
modeling, and will not further deal with RS models; th e 
latter are well covered for example by iMaron 1 (119981) . 
Instead, Sec. [V] will present theoretical case studies with 
model descriptions of a few examples of nanofrictional 
phenomena, such as electronic friction, magnetic dissipa- 
tion, carbon nanotubc tribology. For a closer and fresh 
perception of where the field stands we will also include a 
partial and incomplete list of problems that, in our per- 
spective, still stand out for future theoretical study in 
friction and nanofriction. 



II. SIMPLE NANOFRICTIONAL MODELS 



A. The Prandl-Tomlinson model 



The Prandtl-T omlinson (PT) model (jPrandtll . Il928t 
lTomlinsonl[l929h . which we discuss here in some detail, is 
the most successful and influential MM so far suggested 
for description of nanoscale friction. In particular it ad- 
dresses friction force microscopy (FFM) where friction 
forces are measured by dragging an AFM tip along a 
surface. Qualitative conclusions drawn with this model 
provide guidance to understand friction at the nanoscale, 
that often retain their validity in more advanced models 
and MD simulations. 

PT assumes that a point mass m (mimicking e.g., the 
AFM tip) is dragged over a one-dimensional sinusoidal 
potential representing the interaction between the tip and 
a crystalline substrate. The point tip is pulled by a spring 
of effective elastic constant K, extending between the tip 
position x and the position of the microscope support 
stage, that is driven with a constant velocity v relative 
to the substrate, see Fig. QJ,. Thus the total potential 
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Figure 1 Stick-slip in (a) a cartoon of the PT model; (b) En- 
ergy landscape for a soft spring (low K). The total potential 
(harmonic spring + sinusoidal substrate) exhibits different 
metastable minima, giving rise to the stickslip behaviour; (c) 
a representative experimental friction pattern, for increasing 
load. Lateral force vs position traces demonstrate transitions 
from smooth sliding (top) t o single (middle) and m ostly dou- 
ble slips (bottom). From (|Medvanik et all l2006t ). Similar 
patterns can be generated within the PT model. 



experienced by the tip consists of two parts: (i) the tip- 
substrate interaction, and (ii) the elastic interaction be- 
tween the tip and the support, and can be written as 

U(x,t) = C/ O cos ^— x^j + yfa- vtf (1) 

where Uq is the amplitude, a is the periodicity of the tip- 
substrate potential. Note that in an AFM experiment the 
real "spring constant" mimicked by K in the PT model 
is not only due to the torsional stiffness of the cantilever 
but includes also the contribution from the lateral stiff- 
ness of the contact. There is no attempt in the model 
to describe realistically the energy dissipation into the 
substrate (Joule heat) and all dissipation is described by 
a viscous-like force —rrvyx , where 7 is a damping coeffi- 
cient. The instantaneous lateral friction force measured 
in FFM experiments reads F = —K(x — vt), and the 
kinetic friction Fk is the time average of F. 

The PT model predicts two different modes for the 
tip motion, depending on the dimcnsionless parame- 
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tcr 77 — 4ir 2 U~o/(Ka 2 ), which represents the ratio be- 
tween the stiffnesses of the tip-substrate potential and 
the pulling spring. When 77 < 1 the total potential U (x) 
exhibits only one minimum and the time-dependent slid- 
ing motion is smooth; for 77 > 1 two or more minima 
appear in U(x), and the sliding is discontinuous, charac- 
terized by stick-slip, Fig. [TJd. The value 77 = 1 represents 
the transition from smooth sliding to slips by one lattice 
site (single-slip regime). 

Physically, stick-slip motion corresponds to jumps of 
the tip between successive minima of U(x), due to in- 
stabilities induced by the driving spring (dU/dx = 
0, d 2 U/dx 2 = 0). Close to the inflection point the height 
of the barrier preventing the tip sliding decreases with 
the applied force, as AE ex (const — F) 3 ^ 2 ( Dudko et all 
l2002t iMalonev and Lackst 120061: ISang et all 1200 ll) This 
type of externally induced topological change in a 
free energy landscape is known fold catastrophe, 

and it has been found in many driven systems, in- 
cludin g superconducting quan tum interference devices 
( GareJ . 1 19951: iKurkiiarvi . 19721) . me chanically deformed 
glasse s (IJohnson and Samwerl I2005T) and stretch ed pro- 
teins (jBerSvichTIairEQia EackTeToE [20l3). The 
simulation results obtained for diverse systems show that 
the fold catastrophe scaling is in fact accurate not only 
in the immediate vicinity of the inflection point but over 
reasonably large intervals of loads. 

The possibility of slips of higher multiplicity (multiplc- 
slip regime) occurs fo r larger values of 77 > 4.604 
(|Medvanik et all |2006[ ). It should be noted that this is 



a necessary but not sufficient condition to observe mul- 
tiple slips, since the observed dynamics depends also on 
the damping coefficient 7. In particular, for 77 > 4.604 
one can distin guish b etween the overdamped regime of 
motion, 7 > ^Uo/m 2ir/a, where the tip jumps between 
nearest-neighbor minima of the potential, and the under- 
damped regime, 7 < y^Uo/m 2ir/a, where the tip may 
perform multiple slips over a number of lattice sites and 
even overshoot the lowest well of the potential U(x). In 
that case the minimal spring force reached during stick- 
slip oscillations is negative. 

The elastic instability occurring for 77 > 1 results in 
a nonzero value of the low-velocity kinetic friction that 
is given by the energy drop from the point of instabil- 
i ty to the next minim um of the potential divided by a 
(jHelman et all Il994l) . For 77 < 1 this instability does 
not exist, friction is viscous, and Ff. — > for v — > 0. 
The emergence of static friction can be interpreted as 
the arousal of a saddle-node bifurcation as a f unction of 
77, rea lizing a sort of fold-catastrophe scenario ( Gilmorel . 
119811) . 

Note that in real systems at finite temperature, hys- 
teresis and dissipation must always disappear in the zero- 
speed limit of adiabatic sliding, where stick-slip insta- 
bilities are preempted by thermal fluctuations. This 
regime, sometimes termed thcrmolubricity, is addressed 
in Sect. III.C1 particularly by Eq. fll}. 

In experiment, the effective value of the PT parameter 



77 can be controlled by the variation of the normal load on 
the contact, which changes the potential corrugation Uq 
more than the contact stiffness. FFM experiments at low 
normal loads indeed demonstrated smooth sliding with 
ultralow fr iction, connected to th e absence of elastic in - 
stabilities (jMedvanik et all 120061: ISocoliuc et all . |2004) . 
At higher loads instead, "atomic" stick-slip took place 
with the atomic periodicity of the substrate lattice, while 
increasing load (corresponding to increasing Uq) further 
led to a multiple slip regime as predicted by the PT 
model, see Fig. QJ. 



B. Extensions of the Prandtl-Tomlinson model 

Several generalizations of the original, one-dimensional 
PT model have marked new steps toward understanding 
and implementation of frictional phenomena. These ex- 
tensions included considerations of: 

• two-dimensional structure of surfaces that led to 
the introduction of f ri ctiona l imag in g of interfaces 
dFusco and Fasolinol 12004 1200a iGvalog et all . 
ll995HPrioli et a/.l . I2003T) : 

• thermal fluctuations that allowed to under- 
stand an origin of velocity dependence of fric- 
tion and introduced a new reg ime of fric - 
tion, named "the r molub ri citv" dDudko et al. 



2002t iGnecco et all [20001: iKrvlov et al. 



Sang et all . 120011) : 



Reimann and Evstigneevl 2004Tlliedo et all 2003 : 



2005 



• coupling between normal and lateral motion o f 
the slider (|Rozman et al.l . ll998atlZaloi et aZ.l . ll999l ) 
that led to a new approach to control friction and 
wear by modulating the no rmal load ( Lantz et all 
l2009t [Socoliuc et al. I. l2006h : 

• flexibility of the AFM tip apex that led to a 
predictions of new regimes of motion exh i biting 
complex stick-slip p atterns ( Krvlov et all 120061 : 
iTshiprut et all 120081) . 

Deferring some of these points bearing contact with 
the Frenkel-Kontorova model, we focus first on the effect 
of temperature on friction. 



C. Thermal and velocity effects on nanoscale friction 

The main aspects of thermal effects on fricti o n wer e 
considered in the pioneering work by iPrandtll (|1928[ ). 
Thermal effects can be incorporated into the model ([1} by 
adding a thermal random force f(t) to the conservative 
force between the slider and substrate and the damping 
term —77172;. Then the tip motion is described by the 
following Langevin equation 



rax + 7717a; 



dU(x,t) 
dx 



■/(*) 



(2) 



The random force should satisfy the fluctuation- 
dissipation theorem; as usual, it is chosen with zero mean 
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/(*) 



and (5-corrclated: 



f{t)f{t'))=2m 1 k B T5{t-t'), 



(3) 



where k B denotes the Boltzmann constant and T tem- 
perature. The random forces and the damping term arise 
from interactions with phonons and/or other fast excita- 
tions that are not treated explicitly. 

In the thermal PT model, Eqs. (J2J and ©, beside 
the PT-parameter 77, thermal fluctuations bring out a 
new dimensionless parameter 8 representing the ratio be- 
tween the pulling rate v/a and the characteristic rate 
of thermally activated jumps over the potential barriers, 
uio exp (— Up/knT), where ujq is the attempt frequency 
(|Krvlov et all |2005|) . As a result one should distin- 



guish between two regimes of motion: (i) 5 <C 1, regime 
of very low velocities or high temperatures (typically 
v < 1 nm/s at room temperature), where the tip has 
enough time to jump back and forth across the barrier, 
and (m) S 3> 1, stick-slip regime of motion, where thermal 
fluctuations only occasionally assist the tip to cross the 
barrier before the elastic instability is reached. In these 
two regimes the follo wing expressions for kinetic friction 
have been sugge s ted dDudko etal\ , 120021 : iKrvlov et all 
l2005l : ISang etal\ . l200lh : ~ 



F k (v, T) = a (T) v + O (v 3 ) , 5 < 1, 



(4) 



F k (v t T) = F bT^ln^^J, (5) 
S > 1 and v < BT. 



Here F$ is the athermal (T = 0) low-velocity limit of 
friction, a(T) cx (K/cjq) (U /(k B T)) exp (U /(k B T)) is 
the equilibrium damping felt by the tip that is indepen- 
dent of the ad-hoc damping coefficient 7, and 6, B are 
positive constants which depend on m, K, a, Uq and 
7 but not on v and T. Equation Q, describing the 
slow friction regime called thermolubricity ( Krvlov et all 
120051 ). corresponds to the linear- response regime, while 
Eq. ([5]) has been derived assuming that thermally acti- 
vated depinning still occurs in the vicinity of the ather- 
mal instability point. The velocity and temperature de- 
pendencies of friction force predicted by Eq. ([5]) result 
from the fold catastrophe scaling of the potential bar- 
riers, AE cx (const — F) 3 ^ 2 , which has been discussed 



in the previous section. Furthermore, in between the two 
regimes described by Eqs. (jH) and one should observe 
a logarithmic dependence of Fk on velocity. However, 
it is very difficult to distinguish between [ln(u)] 2 ^ 3 and 
simple In (v) behavior in experiments and numerical sim- 
ulations (|Muser[ [20lTh . The logarithmic (or [ln(u)] 2/3 ) 
regime tends to span many decades, until v becomes so 
large that the inertial or viscous-like effects set in. The 
[ln(w)] 2 ^ 3 dependence of the average rupture force has 
been also found in single-molecule unbinding experiments 



where the energy landscape of complex b iomolecules is 
probe d by applying time-dependent forces ( Dudko et all 
l2003h . 

The theoretical framework outlined above has ex- 
plained a numbe r of FFM experimen t al results on single 
crystal surfaces dGnecco et all l2000l iRiedo et~al. 1. 120031: 
IStills and Overnevl 120031 ). Furthermore, the statistical 
distribution of friction forces w as measured to ma t ch pre - 
dictions from the PT model ( Schirmeisen et al\ , l2005f) . 
These results provide strong evidence that atomic stick- 
slip in FFM is attributable to thermally activated slip 
out of a local minimum as described by the PT model. 
Thermally activated stick-slip friction is only seen in MD 
at sufficiently low speeds, which are s o far on ly achievable 
through accelerated MD (|Li et all l2011al ). At higher 
speeds, friction is mostly determined by dissipative ather- 
mal dynamical processes, which correspond to a funda- 
mentally different regime of sliding. This limits severely 
the regime of validity of comparisons of the PT model 
with MD simulations. 

Equations ^ and also predicts that kinetic 
friction shoul d decr e ase with in c reasin g temperature 
dDudko et all l2002t ISang et all 1200 it ISteiner et all . 
I2009h . Thermal excitations in fact help overcome en- 
ergy barriers and reduce the stick-slip jump magni- 
tude, so that nanofriction should decrease with tem- 
perature provided no other surf ace or material parame - 
ters are altered by temperature ( Szlufarska et aLl . l2008|) . 
Up to now, most FFM measurements have been per- 
formed at room temperature, so that the temperature 
dependence of nanoscale friction has rarely been ad- 
dressed in experim e ntal w or k. Recent experimental 
results (iBarel et al. I, l2010allbl : ISchirmeisen et all 120061: 
IZhao et all l2009t) . however, strongly disagree with the 
predictions of Eqs. (|4]) and (0. Friction forces exhibit 
a peak at cryogenic temperatures for different classes 
of materials, including amorphous, crystalline, and lay- 
ered surfaces. Can this effect be explained within the 
PT model? Recent analysis of the thermal P T model 
(jFaiardo and Mazol [20lbTlTshiprut et ^.1 . 120091 ) demon- 



strated that the friction force may indeed exhibit a peak 
in the interval of temperatures corresponding to a tran- 
sition from a multiple-slip regime of motion, at low T, 
to the single-slip regime at higher T. In this picture, 
interplay between thermally activated jumps over poten- 
tial barriers and the reduction of the slip spatial exten- 
sion with T may lead to a nonmonotonic temperature 
dependence of friction. However, the PT model fails to 
reproduce the observed features of the temperature and 
velocity dependencies of kinetic friction, and of the force 
traces measured with atomic resolution. 



D. The Frenkel-Kontorova model 

The basic model describing the sliding of crystalline in- 
terfaces is t he one-dimensional F r enkel -Kontorova (FK) 
model, see iBraun and Kivsharl <j2004h and references 
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Figure 2 A sketch of the FK model with the two competing 
lengths: interparticle and substrate periodicities. 



therein. First analytically treated by iDehlingerl ( 19291 ) 
and then introduced to descr i be dislocations in solids 
(iFrenkel and Kontoroval Il938t iKontorova and Frenkell . 
lL938al lbh. the FK model found subsequently a wide area 
of applications, in particular, in surface physics, where it 
is often used to unravel the physical behavior of adsorbed 
monolayers, specifically to address competing incommen- 
surate periodicities. 

The standard FK model Hamiltonian is 



Pi K 

2^ + T (x * +1 



U 2-KXi 

— cos 

2 a b 



,(6) 



describing a ID chain of N harmonically coupled classical 
"atoms" subjected to a sinusoidal potential, see Fig. [3] 
The first term in Eq. © is the kinetic energy of the 
chain, the second one describes the harmonic interaction 
of the nearest neighbors in the chain with elastic constant 
K and equilibrium distance a c , and the last term is the 
interaction of the chain particles with the periodic poten- 
tial of magnitude Uq and periodicity a&. Static friction is 
probed by driving all atoms with an extra adiabatically 
increasing force F until sliding initiates. 

Tribological processes in the FK model are ruled by 
kink (topological soliton) excitations. Consider the sim- 
plest case of the trivial commensurate ground state when 
the number of atoms N coincides with the number of 
minima of the substrate potential M, so that the dimen- 
sionless concentration 8 = N/M = at/a c is 1. In this 
case, adding (or subtracting) one extra atom results in 
a chain configuration with a kink (or an antikink) exci- 
tation. After relaxation, the minimum-energy configu- 
ration corresponds to a local compression (or extension 
in the antikink case) of the chain. Kinks arc important 
because they move along the chain far more easily than 
atoms: the activation energy for kink motion (the Peierls- 
Nabarro barrier) is always smaller or much smaller than 
the amplitude Uq of the substrate potential. The rel- 
evance of topological defects for slip can be understood 
going back to the pioneering work by Frenke l on the shear 
strength of crystalline solids ([Frenkell . Il926|) . Frenkel es- 
timated the ideal plastic yield stress of a crystal as the 
stress needed to displace one atomic plane by one lat- 
tice spacing. This calculation leads to an estimate for 
the yield stress which is of the order of the shear mod- 
ulus. The result is in contradiction with experiments 
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Figure 3 Detailed behavior (atomic trajectories versus time) 
at the depinning transition at a small nonzero temperature of 
the FK chain with 9=1. The onset of motion is marked by 
the creation of one kink-antikink pair. The kink and antikink 
move in opposite directions, collide quasielastically (because 
of the periodic boundary conditions), and soon a second kink- 
antikink pair is created in the tail of the primary kink. This 
process repeats with an exponential (avalanche-like) growth of 
the kink-antikink co ncentration, leading to the totally sliding 
state. Adapted from iBraun et al\ i|l997f ). 



which typically show much smaller results. The reason 
for this discrepancy is rooted in the presence of disloca- 
tions which can be displaced by much smaller stresses, 
leading to plastic deformation much earlier than in the 
perfect-crystal limit. In ideal conditions, the only bar- 
rier to dislocation motion is the Peierls-Nabarro stress 
due to the periodic lattice and this is typically orders of 
magnitude smaller than the shear modulus. 

Because the kinks (antikinks) correspond to extra 
atoms (vacancies), their motion provides a mechanism 
for mass transport along the chain and are thus responsi- 
ble for mobility, conductivity, and diffusivity. The higher 
the conce ntration of kinks, the higher will be the system 
mobility ( Vanossi et all , |2003|) . When the ground state 
is commensurate (i.e., 6 = 1), at nonzero temperature, 
the first step to initiate motion in the FK model is the 
creation of a kink-antikink pair, see Fig. [3] 

When the elastic layer is of finite extension, kinks are 
usually generated at one end of the chain and then prop- 
agate along the chain until disappearing at the other free 
end. Each run of the kink (antikink) through the chain 
results in the shift of the whole chain by one lattice con- 
stant ab- In the case of a finite film confined between two 
solids, one may similarly expect that the onset of sliding 
is initiated by the creation of a local compression (kink, 
misfit dislocation) at the boundary of the contact, while 
kink's motion is the basic mechanism of sliding. 

A crucial role in the FK model is played by incommen- 
surability and the Aubry transition ( Pevrard and Aubrvl . 
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119831) connected with it. Let the substrate period a b and 
the natural period of the chain a c be such that, in the 
limit of infinite system's length, their ratio 8 = a b Ja c is ir- 
rational. Roughly speaking, in this case the FK chain ac- 
quires a "staircase" deformation, with regions of approx- 
imate commensurability separated by regularly spaced 
kinks (or antikinks if 9 < 1). If there is a nonzero proba- 
bility to find particles arbitrarily close to the maximum 
potential energy Uq these kink s are unpinned and mo - 
bile, otherwise they are pinned ( Floria and Mazol Il996|) . 
For a fixed amplitude of the substrate potential U , the 
FK ground state undergoes a transition between these 
two states (the Aubry transition) at a critical value 
K = K c of the chain stiffness. K c depends dramati- 
cally and discontinuously on the incommensurability ra- 
tio ab/a c defining the interface. In particular, it has been 
proven that K c takes the minimal possible value equal 
to ps 1.0291926 [in units of 2U {n/a b ) 2 } for the ratio 
equal to the irration a l gold en mean a b /a c = ( 1 + \/5) / 2 
(jBraun and Kivsharl . \2004 . From a physical point of 
view, this means that for K > K c there is a continuum 
set of ground states that can be reached adiabatically by 
the chain through nonrigid displacements of its atoms at 
no energy cost (sliding mode). On the other hand for 
K < K c , the atoms are all trapped close to the minima 
of the substrate potential and thus require a finite energy 
per kink (equal to the Peierls-Nabarro barrier) to move 
over the corrugated substrate. Thus, for incommensurate 
contacts above the Aubry transition (K > K c ), the kinks 
are mobile, chain sliding is initiated by even the small- 
est driving force and, accordingly, the static friction force 
vanishes, F s = - the chain sliding is superlubric. On 
the other hand, below K c the two incommensurate ID 
surfaces are locked together due to pinning of the kinks 
that separate local regions of common periodicity, and in 
this case we expect stick-slip. 

The kinetic friction p roperties of the FK model 
( Strunz and Elmerlll998al lbT) are probed by adding a (e.g. 
Langevin) thermostat as described for the PT model 
above. Even where (above the Aubry transition) F s = 
the kinetic friction force Ff. is nonzero, because the dy- 
namics at any finite speed results in the excitation of 
phonons in the chain. Note also that a finite-size T = 
FK model is always statically pinned, even for an ir- 
rational value of a b /a c because of the locking of the 
free ends of the chain. However an Aubry-like transi- 
tion, exhibi ting a symme t ry-breaking nature, can still 
be defined dBenassi et oE 1201 lbt iBraiman et aR Il990t 
iPruttivarasin et all |201ll ). At finite T, pinning can 
be overcome by thermal fluctuations, which can initi- 
ate sliding even in the most-pinned state, the fully com- 
mensurate one, see Fig. [3] Finally, we remark that 
friction in the dynamically driven FK model describes 
fairly just the onse t of sl iding of a crystalline contact 
(jHammerberg et all Il998l) . while it cannot account for 
the highly inelastic plastic or quasi-plastic deformations 
of the surfaces characterizing real-life friction experi- 
ments. 
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Figure 4 The data points show the average friction fo r ce ver - 
sus the rotation angle measured by iDienwiebel et al\ (|2004l ). 
The curve through the data points shows the calculated fric- 
tion force from a generalized PT model for a symmetric 96- 
atom flake. From IVerhoeven et ali (|2004l ). 



E. Superlubricity 

Supcrlubricity is the phenomenon in which two incom- 
mensurate periodic surfaces may slide in dry contact with 
no atomic scale stick-slip instabilities which, as discussed 
above, are the main source for energy dissipation. Its 
physical origin is first of all that the energy of two inter- 
acting infinite incommensurate systems is independent of 
their relative position; and secondly that if they are hard 
enough they will slide without stick-slip. 

Vanishing static friction has been first obtained 
within the FK m o del i n the pioneering work of 
iPevrard and Aubrvl ( 1983f ) for mutually incommen- 
surate periodicities, and sufficiently hard infinite 



lattices. Later , Hirano and Shinid (I1993L Il990l) : 



IShinio and Hiranol (| 19931) predicted that for infinite 



commensurate contacts also the kinetic friction should 
vanish, and called this effect superlubricity. In these con- 
ditions, the lateral corrugation forces between two non- 
matching, rigid crystals cancel out systematically, so that 
the kinetic friction of externally driven solid vanishes at 
zero speed, and is dramatically reduced even at finite 
speed. 

The term superlubricity has been criticized as mislead- 
ing, since it might wrongly suggest zero friction in the 
sliding state in analogy to superconductivity and super- 
fluidity. Instead, incommensurability of periodic inter- 
faces cancels only one of the channels of energy dissipa- 
tion, that originating from the low-speed stick-slip insta- 
bility. Other dissipative processes, such as the emission of 
sound waves, still persist, and therefore even in the case 
of complete incommensurability the net kinetic friction 
force does not vanish. Nonetheless, in the superlubric 
regime one expects a substantial reduction of the friction 
force relative to a similar, but commensurate case. 
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Figure 5 The modified PT model used in the simulations 
of superlubricity. A rigid flake consisting of N atoms (here 
N = 24) is connected by an x spring and a y spring to the 
support of the microscope. The support is moved in the x 
direction. The substrat e is modelled as an infin ite rigid single 
layer of graphite. From IVerhoeven et al\ (|2004l ). 



Detailed experimenta l studies of super l ubrici t y wer e 
recently performed by iDienwiebel et al. I (120051 . I2004T) : 
IVerhoeven et al\ ( 20041 ). who measured friction between 
a graphite flake attached to the FFM tip and an atom- 
ically flat graphite surface. Super-low friction forces 
(< 50 pN) are found for most relative orientations of the 
flake and the substrate, for which the contacting surfaces 
find themselves in incommensurate states (see Fig. @|. 
For narrow ranges of orientation angles corresponding to 
commensurate contacts, stick-slip motion was observed 
and friction was high (typi c ally 2 5 pN). A few earlier 
exper iments ( Hirano et all 1 1 9 9 ll ISheehan and Lieberl . 
Il996l ) also provided indications of superlubricity in dry 
friction. 

These observations of superlubricity can be described 
within a generalized PT model treating the graphite flake 
as a rigid finite lattice, composed of hexagonal carbon 
rings, as shown in Fig.[S] The interaction potential is ob- 
tained by summing the pairwisc interactions between car- 
bon atoms in the flake and in the graphite surface. The 
resulting flake-surface potential £/fl a k G {x c , y c , 4>) depends 
on the position of the center of mass of the rigid flake 
given by the two-dimensional coordinate r c = (x c ,y c ), 
and on the orientational (misfit) angle (j> of the flake rel- 
ative to the surface lattice. The motion of the flake at- 
tached to the FFM tip and driven along the surface is 
described by the PT equation @ where the sinusoidal 
tip-surface potential is replaced by the more complex po- 
tential [/flake (x c , y c , 4>). Assuming that the flake is rota- 
tionally locked (i.e. (j> is constant) the angular dependence 
of average friction force is in agreement with observa- 
tions: the friction exhibits narrow peaks of high friction 
with stick-slip motion around the values of <fi correspond- 
ing to the commensurate configurations of the flake and 
the surface, which are separated by wide angular inter- 
vals with smooth-sliding ultra-low friction corresponding 



to incommensurate configur ations ( Merkle and Marks! . 
120071 ; IVerhoeven et all l2004h . ft should be noted that 
the angular width of the friction maxima, A$, should 
depend on the flake size, tan(A$) = 1/D , where D is 
the flake diameter, expressed in lattice spacings. Accord- 
ingly, the width of friction peaks can be used to estimate 
the flake diameter. 

Superlubricity between incommensurate surfaces pro- 
vides a desired low-friction state essential for the func- 
tion of small-scale machines. However, some experiments 
show that flake superlubricity has a finite lifetime: it 
disappears due to a reorientation of the flake into the 
commensurate state ( Filippov et all 20081) as observed in 
a generalization of t he PT model (|Filippov et all l2008t 
Ide Wiin et ^ll2010l) and in tight-binding atomistic sim- 
ulation (|Bonelli et aZ.I . I2009T) . 

Studies of superlubricity may have important impli- 
cations for understanding the macroscopic properties 
of graphite and other solid la mellar lubricants which 
are common s o lid lubricants (iHeimberg et all l200lt 
iRapoport et all Il997l : ISingerJ . Il998fh Mainly used as 
flaky powder, they are applied where liquid lubricants 
cannot be used, and show remarkable nanotribological 
properties which are still not un derstood. Recent MD 
simulations ( de Wiin et aZ.I . [201lh demonstrated that two 
surfaces lubricated by mobile, rotating graphene flakes 
may exhibit stable superlubric sliding as for ideally in- 
commensurate contacts and for surfaces covered by ran- 
domly oriented pinned graphene patches. Under humid 
conditions, the multi-domain surface structures can form 
spontaneously due to the capillary forces which fix ran- 
domly oriented flakes at the sliding surfaces, while in 
vacuum graphite patches are free to reorient to a high- 
friction and high- wear regime. This may provide an an- 
swer to the long-standing problem of why graphite is such 
a bad lubricant i n vacuum, and needs the humidity of air 
to perform well ( Savagd . Il948f ) . 



F. Extensions of the Frenkel-Kontorova model 

Many relevant generalizations of the FK model have 
been proposed so far to cover a large class of frictional 
relevant phenomena; they mainly consist of modifica- 
tions of model interactions or of dimensionality. For 
realistic physical systems (as, e.g., atoms adsorbed on 
a crystal surface), anharmonicity c an be introduced in 
the c hain interatomic potential, see ( Braun and Kivsharl 
2004). The main novelties here include effects such 
as a broken kink-antikink symmetry, new types of 
dynamical solitons (supersonic waves), a breakup of 
the antikink soliton followed by a chain rupture, a 
changed kink-kink interaction. Likewise, nonsinusoidal 
periodic substrates, c haracterized, e.g., by sharp bot - 
toms and flat barriers ( Pevrard and Remoissenetl . 1 19821 ). 



have been investigated to address atoms adsorbed 
on sim ple metal surfaces. Complex unit cell sub- 
strates ( Remoissenet and PevrardL Il984l : IVanossi et all 
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2003), as we l l as quasiperiodic ( van Erp et all . Il999t 
Vanossi et al. . l2Q00h and disordered corr ugated profiles 



(|Cule and Hwal . ll996l : lGuerra et oZl . l2007t ) have also been 
considered in simulations. These deviations from the 
standard FK potential may lead to qualitatively differ- 
ent excitations such as different types of kinks, phoiion 
branches, changes in kink-antikink collisions. From a tri- 
bological point of view, different types of sliding behav- 
ior are to be expected at low-driving forces, when the 
dynamics is mainly governed by the motion of kink-like 
structures. 

An important and more realistic generalization of the 
standard FK chain with relevant consequences for the 
resulting tribological properties (critical exponents, scal- 
ing of friction force with system size, mechanisms of de- 
pinning, etc.) involves increasing the dimensionality of 
the model. Especially the FK 2D generalized versions 
of the model ( Braun and KivsharL[200l lP crsso nl l2000al ) 
are naturally applicable to the description of a contact 
of two solid surfaces (i.e., the case of "dry" friction), 
in particular as is realized in QCM experiments, where 
2D monoatomic islands of ads orbate atoms slide over a 
periodic crystalline substrate ( Krim and Widoml Il988h . 
These approaches are especially powerful in the investi- 
gation of the transient behavior at the onset (or stopping) 
of sliding, which is q uite difficult to stud y in fully realistic 
3D models, see e.g. (|Braun et oUl200lh . 

Noncontact A FM tips osci l lating on top of kink-like ad- 
sorbate regions ( Maier et all . 120081 ) dissipate significantly 
more than near in-registry regions. This mechanism is 
explaine d by the higher softness and mobility of solitonic 
regions (Bennewitz et all. 120001: I Gauthier and Tsuk ada , 
120001 iHoffmann et aUl200lllLoppacher et a?.l . l2000[) . and 
it has been demonstrated by the dynamics of an incom- 
mensurate FK ch ain, forced and p robed by a locally- 
acting oscillation ( Negri et~al\ , |20 10f) . 

In investigating confined systems und er shea r, FK- 
like models with ju st one particle ( Miiserl 120021 : 
iRozman et al. 



( Braun et all 



1996allbh or an int e racting atom ic chain 
20051: IRozman et all Il997l il998bh embed- 



ded between two competing substrates have led to un- 
cover peculiar tribological phenomena related to stick- 
slip dynamics or to t he appearance of "quantized" slid- 
ing regimes of motion dCastelli et all 2009 M anini et all 



120071: Tsantoro et all 120061: I Vanossi e* a/1 120071 120061) 



While some of these phenomena, such as chaotic and in- 
verted stick-slip motion, two types of smooth sliding and 
transitions between them, have bee n already observed 
(iDrummond and Israelachvilil l200lb iDrummond et all . 
12003b . others are still waiting for experimental confir- 
mation. 

Last but not least, the combined Frenke l -Kontorova - 
Tomlinson (FKT) model (jWeiss and Elmerl Il996l Il997| ) 
has been introduced including harmonic coupling of the 
interacting chain atoms to a sliding body. This approach 
resembles the Burridgc-Knopoff model, to be detailed in 
Sec. IIV.B1 where, however, an on-site interaction with 
the lower body replaces a phenomenological dry friction 
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Figure 6 (color online). Sketch of a typical MD simulation of 
a boundary-lubricated interface under shear. Periodic bound- 
ary conditions are applied in the x — y directions. 



law (usually a velocity- weakening one) . The FKT model 
introduces more degrees of freedom than the PT model, 
and it has been used to describe effects of finite size and 
stiffness of the AFM tip and of normal loa d on friction 
(jlgarashi et all . 120081 : iKim and Fail l2009f ). The latter 
effect has been modeled assuming a linear dependence 
of the amplitude Uq of potential corrugation on the ap- 
plied normal force. The validity o f the FKT model ha s 
been tested by 3D MD simulations (|Kim and FalkLl2009j l. 
which confirmed the outcome of the model for most of in- 
vestigated regimes except the limit of very low stiffness 
and high normal load. Unlike the FKT model in which 
the breakdown of supcrlubricity coincides with the emer- 
gence of the metastable states, in the 3D model some 
metastable states appear to reduce frictional force lead- 
ing to nonmonotonic dependence of force on normal load 
and tip compliance. 

Increasing dimensionality and adding realistic features 
to the FK model brings its extensions into closer and 
closer contact to full-fledged MD simulations. 



III. MOLECULAR-DYNAMICS SIMULATIONS 

The simple low-dimensional MMs discussed in Sect . [TT] 
are useful for a qualitative understanding of many phys- 
ical aspects of friction. To address subtler features, such 
as the temperature dependence of the static friction of a 
specific interface or the Joule-heat dissipation, one should 
go beyond MMs including atomistic structural details of 
the interface. Such an approach is provided by MD sim- 
ulations. 

Advances in computing hardware and methodology 
have dramatically increased our ability to simulate fric- 
tional processes and gather detailed microscopic infor- 
mation for realistic tribological systems. MD simula- 
tions are used extensively in sliding nanofriction to pro- 
vide unique insight into the relevant processes, sometimes 
overturning conventional wisdom. They represent con- 
trolled computational experiments where the dynamics 
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Figure 7 (color online). A simulated truncated-octahedron 
AU459 cluster sliding with one of its (11 1 ) face ts over a mobile 
graphite substrate. From [Guerra et ~al\ ( 20101 ). 



of all atoms is obtained by solving numerically Newton 
(or Langevin) equations of motion based on suitable in- 
terparticle interaction potentials and the corresponding 
interatomic forces. The geometry of the sliding interface 
and the boundary conditions (e.g. as sketched in Figs. [5] 
and [3) can be chosen to explore friction, adhesion, and 
wear. A thermostat, or other form of damping, is in- 
troduced in order to eliminate the Joule heat to obtain a 
frictional steady state. Finally, after specifying the initial 
coordinates and velocities of the particles, the classical 
differential equations of motion are integrated numeri- 
cally. 

A worthwhile guide to atomistic MD simulations of 
frictional processes focusing on fundamental technical 
aspects (realistic construction of the interface, appro- 
priate ways to impose load, shear, and the control of 
temperature) can be foun d in th e review artic les by 
iRobbins and Muserl (|200ll) and by iMuserl (|2006l) . For 
the general classical MP approach, we r e fer th e reader 
to the textbooks by Allen and Tildeslevl ( 199lf ) and by 
iFrenkel and Smitl (11996D . 

By following the Newtonian dynamics of a system ex- 
ecuting sliding for a significant amount of time, quan- 
tities of physical interest such as instantaneous and av- 
erage friction force, mean (centre-of-mass) slider veloc- 
ity, heat flow, and correlation functions are numerically 
evaluated. Unlike standard equilibrium MD simulations 
of bulk systems, frictional modeling inherently involves 
non-equilibrium conditions and a non-linear dissipativc 
response to the external driving. A standard practical 
assumption is to add Langevin terms to Newton's equa- 
tions, like in Eqs. ©, © for the PT model at finite 
temperature. We will return later to this point. 

The choice of the appropriate interaction forces 
between atoms represents a major problem. If 
U{Ri, i?2, ■■■Rn} is the total energy of the system (here 
say, the slider plus the substrate), as a parametric func- 
tion of all atomic coordinates {Ri}, the force on atom 
i is Fi = — Vfl^C/, perfectly determined once U is 
known accurately. Unfortunately this is generally not 
the case, for U is determined by the quantum mechan- 
ics of electrons - a much bigger and unsavory problem 
to solve. Ab-initio MD, e .g. of the Car-Parrinello type 
( Car and Parrinellol Il985l) . has not really been of use 
so far in sliding friction, mainly because it can han- 



dle only rather small systems, typically hundreds of 
atoms, for relatively short times, typically < 1 ns. Most 
MD frictional simulations are therefore based on reason- 
able empirical interatomic forces ("force fields"), rang- 
ing from relatively sophisticated energy surfaces account- 
ing for electrons at the density-func tional level or at the 
tight-binding level dXu et all Il992j ). to angle-dependent 
many-particle potentials, to simple pairwise potentials 
(e.g. Lennard- Jones), to basic simple models of clastic 
springs, extensions of FK-type formulations. In practice, 
several reliable force fields, parameterized to fit differ- 
ent ranges of experimental data a nd material combina 
tions, are available in the litera t ure (IBrenner et al\ . 
Garrison and Srivastaval 19951 : Ghiringhelli et all 



1998; 



2005 



Los et all 120051: IWeiner et all Il986l) . While this allows 
qualitative atomistic simulations of sliding friction, it is 
often far from quantitative. For example, during such 
a violent frictional process as wear, atoms may change 
substantially their coordination, their chemistry, some- 
times their charge. Once a specific system is understood 
after the elaborate development of satisfactory poten- 
tials, the mere change of a single atomic species may 
require a painful reparamctcrization of the interatomic 
forces. As a result, systematic frictional studies are a 
tour de force if no suitable set of consistent potentials is 
already available. A promising approach consists in the 
use of the so-called reactive potentials ( Brenner et all 
120021: Ivan Duin eiail , l200lfc IStuartl . |2000|) , capable of de- 
scribing chemical reactions and interface wear, with the 
advantage, for large-scale atomic simulations, of a good 
computational efficiency compared to first-principle and 
semi-empirical approaches. 



A. Thermostats and Joule heat 

In a tribology experiment, mechanical energy is con- 
verted to Joule heat which is carried away by phonons 
(and electrons in metals). In a small-size simulation, 
the excitations generated at the sliding interface prop- 
agate and crowd-up into an excessively small region of 
"bulk" substrate, where they are back reflected by the 
cell boundaries, rather than properly dispersed away. 
To avoid overheating and in order to attain a frictional 
steady state, the Joule heat must therefore be steadily re- 
moved. If this removal is done by means of standard equi- 
librium thermostats such as velocity rescaling or Nose- 
Hoover or even Langevin dynamics, an unphysical dissi- 
pation is distributed throughout the simulation cell, so 
that simulated atoms do not follow their real conserva- 
tive motion, but rather execute an unrealistic damped 
dynamics whic h turns out to aff e ct the overall tribologi- 
cal properties ( Tomassone et all Il997j) . Similarly in the 
PT and FK models, the damping parameter 7 is known 
to modify kinetic and frictional properties, but there is 
no clear way to chose the value of 7. To a lesser or larger 
extent, this lamentable state of affairs is common to all 
MD frictional simulations. 
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Figure 8 (color online). Sketch of a MD simulation of friction. 
To account properly for heat dissipation, the infinitely-thick 
substrate is divided into three regions: (i) a "live" slab com- 
prising layers whose atomic motion is fully simulated by New- 
ton's equations; (ii) a dissipative boundary zone, coincident 
with the deepmost simulated layer, whose dynamics includes 
effective damping (e.g., non-Markovian Langevin-type) terms; 
(iii) the remaining semi-infinite solid, acting as a heat bath, 
whose degrees of freedom are integrated out. 



To solve this problem, one should attempt to modify 
the equations of motion inside a relatively small simu- 
lation cell so that they reproduce the frictional dynam- 
ics of a much larger system, once the remaining vari- 
ables are integrated out. A recent implementation of a 
non-conserva tive dissipat i on sch eme, based on early for- 
mulatio ns by |MagalinskH (Il959l), and subsequent deriva- 
tions b y Kantorovichl (120081 ); iKantorovich and Rompotisl 
(120081 ); iLi and El (|2007l ). has demonstrated the correct 
disposal of friction-generated pho nons, even in the rel- 
ative ly violent stick-slip regime dBenassi et all l2010l 
l2012h . All atoms near the sliding interface follow plain 
vanilla Newton's equation, while the atoms in the deep- 
most simulated layer, representing the boundary layer in 
contact with the semi- infinite heat bath (see Fig. [SJ, ac- 
quire additional non-conservative (and non-Markovian) 
terms which account for the time history of this layer 
throu gh a memory kernel ( Kantorovichl 120081 ; ILi and El . 
120071 ). Nanofriction simulations exploiting this dis- 
sipative scheme have recently been implemented that 
conceptually and practically improve over a traditional 
Langevin simulation. Improvement is achieved in partic- 
ular by adjusting the damping 7 applied to the simulation 
boundary layer so as to variationa lly minimize the ener gy 
back-reflected by that boundary ( Benassi et aZj . [2012h . 



B. Size- and time-scale issues 

Modern CPUs perform of the order of 10 9 floating- 
point operations per second (FLOPS) per core. Classical 



MD can take advantage of medium-scale parallclization, 
with fairly linear scaling to approximately 10 2 cores, thus 
affording about 10 11 FLOPS. As the calculation of the 
force acting on each atom (usually the dominating step in 
a MD calculation) requires, depending on the complexity 
and range of the force field, about 10 to 10 2 FLOPS, the 
product of the number of simulated particles N times 
the number of time-integration steps N ate p per runtime 
second on a modern medium-size parallel computer is 
approximately N N stC p — 10 10 . With a typical time-step 
in the fs range, a medium-size simulation involving N = 
10 5 atoms can progress at a rate of 10 5 fs per second, i.e. 
approximately 10 9 fs = 1 /its in a simulation day. This 
total time scales down for larger systems sizes. 

These estimates should be compared with the typical 
times, sizes, and speeds of tribology experiments. If we 
wish to address macroscopic sliding experiments, typical 
speeds would be in the 0.1 to 10 m/s range: in 1 /its the 
slider advances by 0.1 to 10 /im, i.e. approximately 10 3 to 
10 4 typical lattice spacings, enough for a fair statistics of 
atomic-scale events (but hardly sufficient to gather signif- 
icant data about phenomena such as the diffusion of addi- 
tives or of wear particles within the lubricant, or step- or 
defect-related phenomena). In a nanoscale FFM experi- 
ment, however, the tip advances at a far smaller average 
speed (i.e. ~ 1 /im/s) and we can simulate a miserable 
~ 1 pm advancement in a typical run, far too short to ob- 
serve even a single atomic-scale event, let alone reaching 
a steady state. Therefore, whenever long-distance cor- 
relations and/or slow diffusive phenomena and/or long 
equilibration times are expected, MMs will perform bet- 
ter than fully atomistic MD simulations. There is nev- 
ertheless so much physical insight to be extracted from 
MD simulations that it makes sense to run them even at 
larger speeds than in AFM of SFA experiments; and in- 
deed, the sliding speed adopted in most current atomistic 
MD frictional simulations is in the 0.1 to 10 m/s range. 

One of the challenging problems for MD simulations 
is to account for the transition from stick-slip to steady 
sliding. In SFA and AFM experiments, stick-slip with 
its associated hysteresis and large friction generally dis- 
appears for speeds larger than ~ 1/im/s, while in MD 
simulations the transition takes place in the m/s range. 
This major discrepancy (up to ~ 6 orders of magni- 
tude in speed!) between simulations and exper i ments , 
has been discussed e.g. b y Braun and Naumovetd (I2006T ); 
iBraun and Roderl ((20021 ) : iLuan and Robbinsl (|2004l ). and 
relates to the effective spring-force constants and mass 
distributions, that are hugely different in the two cases, 
and much oversimplified in simulations. Several attempts 
to fill these gaps rely on methods, including hyperdynam- 
ics, parallel-replica dynamics, temperature-accelerated 
dynamics, and on-the-fly kinetic Monte Carlo devised 
in recent years faim and Falkl , 1201 it iMishin et al. I. I2007t 
IVoter et q/.l . l2002D . 

Another important aspect present in experiments and 
largely missed by MD simulations is the ageing of con- 
tacts due to the substrate relaxation. Ageing can de- 
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crease substantially the critical velocity for the tran- 
sition from stick-slip to steady sliding. Contact age- 
ing is also believed to be responsible for the increase 
of the static friction force as a function of the contact 
time. Direct imaging of contact regions in samples un- 
der a normal load show a log arithmic growth with time 
( Dieterich and Kilgorel Il994l) . leading therefore to an in- 
creasing static friction. At the phenomcnological level, 
frictional ageing is well described b y rate and st ate fric- 
tion laws, widely used in geophysics ( Ruinal . [l983|) . but its 
microscopic origin is still debated. The main mechanisms 
that h ave been invoked in t he past to explain it are plastic 



creep (IHeslot et all 119941) or chemical strengthening at 



the interface ( Li et all 1201 lbl ). In a recent paper, AFM 
was used to explore ageing in nanoscale contact inter- 
faces, finding supporting evidence for the second mech- 
anism, since when the c ontact surface w as passivated it 
showed no more ageing (|Li et aLl . [20lTbh . It is however 
likely that at larger scales and loads plastic creep would 
also play an important role. Beyond its direct relevance 
for friction, the intriguing issue of contact ageing occurs 
in other non-equilibrium disordered systems such as gran- 
ular media or glasses. 



Figure 9 (color online). A snapshot of a MD simulation of 
the 2D solitonic pattern in the boundary layer of a solid lubri- 
cant (clear grey) in contact with a perfect crystalline surface 
(dark/red), induced by a 16% lattice-constant mismatch. The 
Lennard- Jones interaction of this simulation favors in-registry 
hollow sites, while unstable top sites mark solitonic regions. 
Other layers were omitted for clarity. 



D. Selected applications of MD in nanotribology 



C. Multiscale models 

Since it is currently impossible to treat atomistically 
all the characteristic length scales that mark the dynam- 
ical processes entering the friction coefficient of engineer- 
ing materials, a rising effort is nowadays devoted to de- 
velop multiscale approaches. The basic consideration is 
that unless conditions are very special, all processes far 
away from the sliding interface can be described at least 
approximately by continuum mechanics, and simulated 
using finite elements, allowing for a macroscopic descrip- 
tion of elastic and plastic deformation. The advantage 
of these continuum-theory methods is that it is possi- 
ble to increasingly coarse-grain the system as one moves 
away from the sliding contact, thereby hig hly reducin; 
the computational effort . Sev eral groups (jLuan et 
120061: iMcGee and Smith! . |2007( ) combine the atomistic 
treatment of the interfacial mating region, where dis- 
placements occur on an atomic or larger length scale, 
with a coarse-grained or finite-element continuum de- 
scription elsewhere, where strains are small and continu- 
ous. The main difficulty lies in the appropriate matching 
conditions betw een the atomistic and continuum regions 
(IE et all 120091) . Since the atomic detail of lattice vi- 
brations (the phonons), which are an intrinsic part of 
the atomistic model, cannot be fully represented at the 
continuum level, conditions must be met, for example 
that at least the acoustic phonons should not be exces- 
sively reflected at the atomistic-continuum interface. In 
other words matching at this interface must be such that 
long wavelength deformations should transmit with rea- 
sonable accuracy in both directions, which is vital to a 
proper account of Joule-heat disposal into the bulk. 



We survey here some recent results from the growing 
simulation literature, mostly from our groups, and cer- 
tainly not providing an adequate review of the field. 



1. Boundary lubrication in confined systems 

Hydrodynamics and elasto-hydrodynamics have been 
very successful in describing lubrication by microscopi- 
cally thin films. With two sliding surfaces well separated 
by a hydrodynamically fluid film, friction is mainly de- 
termined by the lubricant viscosity. The friction coef- 
ficient can be calculated using the Navier-Stokes equa- 
tions, showing a monotonic increase with the r elativ e 
sliding velocity between the two surfaces ( Szeril . l200ll) . 
For small driving velocity and/or high load, the lubri- 
cant cannot usually keep the surfaces apart and solid- 
solid contact eventually ensues. Even before full squeeze- 
out, a liquid confined within a nanometer-scale gap ceases 
to behave as a structureless fluid. Pioneering studies of 
confined systems under shear reveal a sequence of dras- 
tic changes in the static and dynamic properties of fluid 
films in this "bounda ry-lubrication" regime. 

SFA experiments ( Yosbizawa and Israelachvilil Il993l ) 
and MD simulations ( Gao et all 1997al lbT) have both 
shown clear upward frictional jumps, in correspondence 
to squeczout transitions from N to N-l lubricant lay- 
ers. The lubricant squeezout for increasing load becomes 
harder and harder, corresponding t o a (near) crystalliza- 
tion of the initially fluid lubricant dPersson and Mugeld . 



l2004lPersson and Tosattil . ll994l:lTartaglino et a/J . l2006[) 
Owing to layer- by- layer squeezout, there will be frictional 
jumps for increasing load. Friction would not necessar- 
ily always jump upward, because restructuring of the 
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solidified trapped lubricant film, and/or a switching of 
the maximum shear gradient from within the lubricant 
layer (possibly accompanied by local melting), to the 
lubricant-slider interface may take place. It could jump 
downward in particular, if lattice mismatch between the 
compressed boundary lubricant layer and the rigid sub- 
strates jumped, say, from commensurate to incommen- 
surate, the latter supcrlubric with a mobile soliton pat- 
tern, as sketched in Fig. [9] Possibilities of this kind are 
presently the object of investigation. 

Simulations also show that the presence of impurities 
or defects between two incommensurate stiff sliding sur- 
faces can, even in a relatively small concentration, lead to 
pinni ng and nonzero static friction ( Miiser and Robbinsl . 
2000). Defects can destroy superlubricity by introduc- 
ing mechanical instabilities, which may occur depending 
on the system dimensionality, the structure and relative 
orientation of the confining walls, and the detail of the 
lubricant-wall interactions. 

MD investigations of a melting-freezing mechanism 
in the stick-slip phenomen ology of boundary-lubricate d 
films were carried out b y Stevens and Robbinsl (|1993|) : 
iThompson and Robbinsl (|l990t) . Various realistic mod- 
els for lubrication layers in very specific contexts 
have been investigated with e x tensive MD simu l ations 
(iBraun and Naumoyetsl. 120061: IChandross et all 120081 
12004 lLorenz et aZ.L l2010allbh . 



2. Sliding of absorbed monolayers 

An ideal experimental setup to investigate the molecu- 
lar origin of friction is provided by the QCM, where atom- 
ically thin systems, usually just one monolayer or less of 
rare-gas is deposited on a substrate resting on a quartz 
crystal. When the crystal surface oscillates strongly 
enough so as to dislodge the adsorbate islands under their 
own inertial force, the sliding friction is revealed by me- 
chanical damping. Depending on the substrates and on 
the system interactions at play, the equilibrium proper- 
ties of these systems may exhibit distinct structural ge- 
ometries, of interest in the field of 2D phase transitions; 
different dynamical phases may show up for the driven 
over-layer, such as pinned solid, slid i ng solid, and liq- 
uid p hases dBruschi et all 120021: iKriml . 120071: iKrim et all 

[nil. 

These experimentally well-characterized systems 
have also been studied theoretically and numerically 
by m olecular dyn a mics s imulations (|Braun et al 
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Cieplak et all . Il994t iGranato and Yinel. 12004 ; 
Persson and Nitzan . 199£ ; ISokoloff and Tomassone , 
19981 ). When the interface is commensurate the 



static friction is too large to allow for a massive slip 
( Cieplak et aUfl994l) : yet slip is observed experimentally 
for a Xe monolayer on a Cu(lll) su bstrate, a system 
that forms a commensurate interface ( Coffey and Kriml 
|2QQ5f ). A possible explanation is a slip which does not 
occur coherently but by a nuclcation process in which 



a bubble slips forward creating a ne w commensurate 
doma in as revealed by MD simulations ( Reguzzoni et all 
|2010D . It is possible to estimate the critical radius of this 
domain and the energy barrier associated with the nucle- 
ation process by the conventional theory of nucleation, 
estimating the value of the relevan t parameters from a 
solut ion of an effective FK model (IBraun and Kivsharl 
l2004f ). Thanks to thermally activated nucleation, with 
the additional help of impurities and disorder acting as 
seeds, the monolayer can slip under lateral forces that 
are much smaller than those expected for a rigid layer 
pinned by the commensurate interface, providing an 
explanation for the experimental results. 

When the adsorbate is incommensurate and hard, soli- 
tons exist already in the ground state, and their free mo- 
tion should in principle permit supcrlubric sliding. If 
however the adsorbate forms islands, as should generally 
be the case, perfect superlubricity is still broken by the 
island's edges, which present a barrier to the entering and 
exiting of solitons, necessary for the island to slide. This 
island edge pinning is currently under study by N. Varini 
and collaborators. 



3. Extreme temperature/speed conditions 

An advantage of MD is that it can address extreme 
or otherwise unusual frictional situations, still unex- 
plored because they may be experimentally difficult to 
realize. One such example are the high "flash tem- 
perature" regimes caused by local Joule heating due 
to wear and other mach ining or braking conditions 
(Bo wden and Taborl . Il950l) . liven in equilibrium but 
at temperatures close to the melting point, the outer- 



most layers of a solid substrate gener ally undergo "sur- 
face melting" (jTartaglino et all 120051) . In these condi- 



tions, AFM nanofriction cannot generally be experimen- 
tally accessed, because the nearly liquefied surface layers 
jump to contact and wet the tip long before it reaches 
nominal contact ( Kuipers and Frenkenl . Il993). However, 
some solid surfaces, such as NaCl(100), do not melt 
(IZvkova-Timan et al. , 2005), thus making for an inter- 
esting, albeit purely theoretical so far, case study. MD 
simulations predict that high-temperature nanofriction 
over such a nonmelting surface would behave very differ- 
ently depending whether the tip-surface contact is "hard" 
or "gentle" . In the first case, the tip plows the substrate 
with wear. The friction coefficient, very large at low tem- 
perature, drops close to the melting point, when the tip 
itself provokes local melting, and moves accompanied by 
its own tiny liquid droplet, precisely as in ice skating. 
For gentle, low-load, wear free sliding on a hard surface 
the opposite is predicted. Here friction, initially very 
small, is expected to increase as temperature is raised 
close to the melting point where the surface, still solid, 
becomes softer and softer due to increasing anharmonic- 
ity, w ith an analogy to flux lines in type II superconduc- 
tors. (jZykova-Timan et aLl . l2007l ). 
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A second example is high-speed nanofriction, as could 
be expected by a tip or a surface-deposited nanoclus- 
tcr (Fig. O sliding at large speed over a smooth crystal 
surface. A speed in excess of 1 m/s is many orders of 
magnitude higher than that of ordinary AFMs or other 
nanofrictional systems, and is as yet unexplored. MD 
simulations, carried out for the test case of gold clus- 
ters sliding over graphite surfaces, show, besides a stan- 
dard low-speed drift sliding regime, the emergence of a 
novel "ballisti c " slid ing regime, typically above 10 m/s 
( Guerra et all l2010f) . The temperature dependence of 
the cluster slip distance and time, measuring its sliding 
friction, is predicted to be opposite in these two regimes, 
high-speed ballistic sliding and low-speed drift sliding. 
The interplay of rotations and translations is crucial to 
both regimes. Simulations show that the two are cor- 
related in slow drift but anti-correlated in fast ballistic 
sliding. Despite the deep difference with drift, the speed 
dependence of ballistic friction is, like drift, viscous, a 
useful result whose validity was not discounted in princi- 
ple, and which it would be interesting to pursue and test 
experimentally 



4. Nanomanipulation: pinning vs diffusion 

AFM manipulation of surface-deposited clusters can 
serve as a useful method to measure the interfacial fric- 
tion of structurally well-defined contacts of arbitrary size 
and material combinations. Here, MD simulations are 
extremely useful in understanding dcpinning, diffusion, 
and frictional mechanisms of clusters on surfaces. In- 
deed, one of the remarkable experimental observations of 
the last decade concerns the unexpected ability of rela- 
tively large metal clusters to execute friction-free motions 
and even long jumps with size and shape conservation 
dBardo tti et al.. 1996 HBrndiar et al. L I2011I: iDietzel et all 
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Figure 10 The effect of oscillations on the lateral force de- 
tected by AFM scanning forward (solid curve) and backward 
(dotted curve) on an atomically flat NaCl surface. An average 
normal load Fn = 2.73 nN was kept constant (a) without a 
bias voltage between the cantilever and the sample holder and 
(b) with an applied AC vol tage with frequ e ncy / = 56.7 kHz 
and amplitude 1.5 V. From ISocoliuc eta/1 l|2006h . 



5. Simulated frictional control 

Exploring novel routes to achieve friction control by 
external physical means is an important goal currently 
pursued in nanotribology. Two methods have recently 
been suggested by simulation: mechanical oscillations 
and phase transitions. 

Mechanical oscillations. Natural or artificially induced 
oscillations obtained by small normal or lateral mechan- 
ical vibrations may, when applied at suitable frequency 
and amplitude ranges, help driving a contacting inter- 
face out of its potential-energy minima, increasing con- 
siderably surface mobility and diffusion, and reducing 
friction and wear. Flexibility and accessibility are the 
main relevant features of this approach, since frictional 
properties can be tuned continuously by the frequency 
and the amplitude of the applied vibrations. This ef- 
fect has be e n demonstrated experimentally with AFM 
(IJeon et all 120061: lLantz et all 120091 iRiedo et all 120031: 



ISocoliuc et al. . 120061). see Fig. [TU1 and in sheared con- 
fined system (|Bureau et al. , 2000; Cochard et al. 1. 120031: 
iHeuberger et al. , 19981) and numerically wi th atomistic 



MD dCapqzza et all 120091: iGao et a71.ll998l) or MM ap- 



prising typically hundreds of atoms, have been repeat- 
edly observed to diffuse on highly oriented pyrolytic 
graphite (HOPG) surfaces with surprisingly large ther- 
mally activated diffusion coefficients already at room 
temperature; a similar behavior was reported also for 
larger antimony clusters. MD simulations of the diffu- 
sive regime have shown the possible coexistence of stick- 
ing periods, a nd of long jumps, reminiscent of so - called 
Levy flights dGuerra et all \20ld, iLewis et all \200& 
iLuedtke and Landman . 119991: iMaruvamal |200| ). The 
sticking lasts so long as the cluster-substrate surfaces are 
orientationally aligned, and the long sliding jumps occur 
when a therma l fluctuation rotates the cluster destroying 
the alignment ( Guerra et q.Z.1 . [2oToh . Further understand- 
ing of the sliding of these dep osited nano-objects will be 
of co nsiderable future value ( Schirmeisen and Schwar3. 
I2009D . 



120081: iPaolicelli et all [20081 12009D . Gold clusters, com- proaches dRozman et aTTll998at iTshiprut et all 120051: 



IZaloi et ~al. , 1999[) . On a larger scale, it has also been 



reported that in sheared granular media experiments 
the stick-slip behavior is also s ignificantly per t urbed 
by tiny transverse vibrations (ICapozza et al\. 2011 
Johnson and Ji al. 120051: [Johnson et all 120081: iKrim et al 



201l[ ). Despite these promising numerical and experi- 



mental contributions, a detailed analysis accounting for 
the friction dependence on vibrations is still to some ex- 
tent lacking. Most past theoretical studies of mechanical 
control adopted an oversimplified single-asperity model, 
which misses the collective behavior of multi-asperity 
mesoscopic interfaces. This and other aspects of inter- 
face oscillation are still calling for a proper treatment 
and understanding. 

Phase transitions. Another idea to control friction is to 
employ a substrate undergoing a phase transition. While 
it is obvious that friction will change in the presence of a 
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phase transition, it is more subtle to qualify and quantify 
precisely the effect. Surprisingly perhaps for such a basic 
concept, there are essentially no experimental data avail- 
able - and no theory either. W hile it would be temptin g 
to use linear-response theory (|Ala-Nissila et all Il992t) . 
with the fluctuation-dissipation theorem making direct 
contact between criticality and viscous friction, one can- 
not ignore that realistic dry friction is dominated by 
stick-slip instabilities that are intrinsically violent and 
non-linear. Hence, one is left with MD simulations. A 
PT-like MD nanofrictional simulation based on a point 
slider over a 2D model substrate with a built-in struc- 
tural displacive transition recently predicted that stick- 
slip friction should actually peak near th e substrate crit- 
ical temperature ( Benassi et all l2011al ). These results 
show that friction will also depend upon the order param- 
eter direction, a dependence due to the different ability 
of the slider to elicit disorder in the substrate depending 
on that direction. Some level of control of atomic-scale 
friction can thus be anticipated through external switch- 
ing of the order parameter direction (e.g., by an external 
field or strain). Although the magnitude of the phase 
transition effects relative to the background friction will 
of course depend on the real system chosen, these re- 
sults suggest pursuing this idea experimentally in, e.g., 
displacive ferro and antiferrodistortive materials, a vast 
realm of solids exhibiting continuous or nearly continuous 
structural transitions. 



IV. MULTICONTACT MODELS 

The PT model discussed in earlier chapters provides 
a good initial description of the frictional behavior of an 
individual contact that can be relevant for nanotribol- 
ogy e xperiments. However, recent simulations (|Mo et all 
I2009D revealed that, even for an apparently sharp AFM 
tip sliding on a crystalline surface, the actual interface 
consists of an ensemble of individual atomic contacts (see 
Fig. ITT]) . On larger scales such as the mesoscale of SFA 
experiments the multicontact picture becomes even more 
obvious. The way how individual contacts can be aver- 
aged to yield macroscopic friction law has been the fo- 
cus of intense research in the past decades. Friction is 
not simply the sum of single- asperity responses, but is 
influenced by temporal and spatial dynamics across the 
entire ensemble of contacts that form the frictional in- 
terface. Long-range elastic interactions b etween contacts 
are i mportant and canno t be neglected. iPersso 3 (l200lh 
and iPersson et al\ { 2005| ) greatly improved and gener- 
alized these concepts to more realistically fractal rough 
surfaces. 




Bond rupturing 



Figure 11 (color online). Multicontact modeling, (a) Sketch 
of a typical geometry for an amorphous tip sliding on a flat 
crystalline surface, (b) A model to simulate multiple contacts 
at the tip-sample interface. Rates of contact formation and 
rupturing processes are determined by the heights of the cor- 
responding energy barriers, AiJ on and A_E g . Adapted from 
iBarel et al\ (|2010al V 



kinetic models (iBarel et all l2010allbl; iBraun et all 12009 ; 



A. Mechano-kinetic models 

A significant progress in the solution of these problems 
has been recently achieved with coarse-grained mcchano- 



Braun and Roderl , 120021: iFilippov et all |2004 IPersson , 
19951 ) that describe friction through dynamical rup- 
ture and formation of interfacial contacts (junctions). 
These contacts may represent molecular bonds, capillary 
bridges, asperities between rough surfaces, and for lu- 
bricated friction they can mimic patches of solidified lu- 
bricant or its domains. Each contact is modeled as an 
clastic spring connecting the slider and the underlying 
surface. As long as a contact is intact (unbroken), it is 
increasingly stretched with a speed equaling the veloc- 
ity of the slider, and thus produces an increasing force 
that inhibits the motion; after the instability point is 
reached, a ruptured contact relaxes rapidly to its mi- 
stretched equilibrium state. The kinetics of contact for- 
mation and rupturing processes depends on the physi- 
cal nature of contacts. For atomic scale contacts, capil- 
lary bridges and domains of solidified lubricants, the pro- 
cesses of rupture and formation of contacts are thermally 
activated, and the interplay between them may lead to 
a complex dependen ce of friction on slid er velocity and 
sample temperature (iBarel et al. L 2010allbl: Braun et al. , 
' 20091; IBraun and Roderl. l2002t IBraun and Tosattil I200S ; 
Filippov et all 12004 iPerssonl ll995l T~ For micro- and 
macroscopic asperities between rough surfaces thermal 
effects are less significant and threshold rupture forces 
should be taken from a distribution that is determined by 
the structure of the contacting surfaces. The mechanism 
of contacts detachment is similar to t he one proposed pr e- 
viously by the fiber-bundle models ( Alava et aLl . l2006h . 
At the nanoscale the rates of formation k on and rup- 
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turing k a g of microscopic contacts arc defined by the cor- 
responding energy barriers AE on and AE a g. The barrier 
for detachment, AE D g, is force dependent and diminishes 
as the force acting on the contact increases and the con- 
tact is stretched. As we discussed above, precisely this 
mechanism characterizes the PT model, but what has 
rarely been rationalized so far is that the process of con- 
tact formation must be considered as well. 

The dynamics of friction in the mechano-kinctic mod- 
els is determined by four characteristic rates: (i) the rate 
k° e of spontaneous detachment of contacts, (ii) the rate 
k on of contact formation, (iii) the rate Kv/f s of forced 
unbinding, and (iv) a characteristic rate of the pulling 

force relaxation, cj m = max ^ Kj (7m), yj K/rnj . Here 

f s is an average threshold force for the contact rupture, 
K is the stiffness of the pulling spring, 7 is the dissipation 
constant, and m is the mass of the slider. Correspond- 
ingly, these models exhibit three regimes of motion: (i) 
smooth sliding at very low velocities or high tempera- 
tures, with fc° ff > Kv/f s , (ii) smooth sliding also for 
high velocities or low temperatures, with fc on < Kv/f s , 
and (iii) stick-slip oscillations for intermediate velocities 
and temperatures. These multicontact models demon- 
strate that the overall smooth sliding corresponds to un- 
correlated atomic-scale stick-slip (or smooth) motion of 
individual junctions, while the global stick-slip motion 
emerges from a cooperative behavior of the junction sub- 
systems. It is interesting to note that the transition from 
smooth sliding to stick-slip with increasing v was indeed 
observed in SFA experiments with two weakly adherin g 
boundary- lubricated surfaces ( Drummond et all 120031 ) . 

An important advantage of the mechano-kinetic mod- 
els is that they are directly scalable to meso and macro- 
scales. Application of these models already allowed to 
solve some significant disagreements between the experi- 
mental observations and the predictions of the PT model 
and MD simulations. Firstly, SFA experiments found 
that the critical velocity for transition from stick-slip is 
in the interval of 1 — 10 /xm/s, while the MD simula- 
tions and PT model lead to values which are 6 o r 7 or- 
ders of magnitude larger dBraun and Naumovetsl. 120061 
iBraun and Rodeil l2002t iLuan and Robbinsl . |2004| ). Ac- 
cording to the mechano-kinetic models the transition 
should occur at v w f s k on /K that for reasonable val- 
ues of the parameters agrees with the experimentally 
observed values of the critical v elocity. Second l y, the 
PT model and MP simulations (iBrukman et all 120081: 
iDong et all 1201 it ISteiner et all , |2009|) fail to reproduce 
the nonmonotonic temperature dependence of the aver- 
age friction force found in FFM experiments for several 
material classes. The mechano-kinetic model demon- 
strates that the peak in the temperature dependence of 
friction emerges from two competing processes acting at 
the interface: the thermally activated format ion and rup- 
turing o f an ensemble of atomic contacts ( Barel et all 
l2010aM . This observation also provides a direct link be- 
tween the temperature and velocity dependencies of fric- 
tion and it shows the experimentally observed fingerprint 



in the friction-velocity curves. Specifically, at tempera- 
tures above the peak temperature, friction increases with 
scan speed, whereas, below the peak, friction decreases 
with velocity. 

An important and still unresolved question is what is 
the microscopic origin of the rupture and reattachment 
processes i ntroduced fo r the in te rpretation of FFM ex - 
periments (|Barel et all l2010allbl : iFilippov et 'ail . 12004) . 
Simila rly to the mechanism of energy dissipation in 
AFM (jGhasemi et aZll2008tlKantorovich and Trevethanl . 
Hqq3), they can be attributed to reversible jumps of sur- 
face atoms, flips of surface fragments or transitions be- 
tween different tip structures, which are induced by the 
tip motion along the surface. These dissipative processes 
result in a bistable potential-energy profile for the tip- 
surface junction where the barrier separating the poten- 
tial minima is continuously changed during sliding. An 
unambiguous understanding of the nature of the corre- 
sponding instabilities and evaluation of microscopic pa- 
rameters which determine the values of the rupture and 
reattachment rates require first-principles calculations of 
potential energy surfaces for the tip-surface junctions. 
The first attempts in bridging the gap between mechano- 
kinctic models and MD simulations have already demon- 
strated that model parameters can be completely speci- 
fied usin g information obta ined from fully atomistic sim- 
ulations ( Perez et a?.l . l2010l ). Then the parameter- free ki- 
netic models are able to reproduce the temperature and 
velocity variation in the friction force as obtained from 
fully dynamical atomistic simulations with very high ac- 
curacy over a wide range of conditions. This combined 
approach is promising because it allows the full atom- 
istic details provided in MD simulations to be used in 
interpreting experimental phenomena at time and length- 
scales relevant to tribological measurements. 



B. Elastic interactions and collective effects 

Elastic instabilities play, as we saw, a prominent role in 
explaining frictional dissipation at the nanoscale and one 
may thus ask what happens at larger scales for multicon- 
tact interfaces. The role of elasticity in friction crucially 
depends on whether stress gradient arc present or not 
( Lorenz and Perssonl 12012]) . For uniform loading, stress 
is distributed homogeneously on the contact surface and 
elastic interactions mediate the response of local slip fluc- 
tuations. Experimental conditions, however, often lead 
to shear stress concentration at the edge of the sample 
leading to detachment fronts. We discuss the first case 
in this section and the second in the the next one. 

The interplay between small scale disorder due to ran- 
dom contacts and elastic interactions between the con- 
tacts is a complex statistical mechanics problem that is 
encountered in different contexts from vortex li nes in su- 
perconductors ( Larkin and Ovchinnikovl 1 19791) . disloca- 



tions in solids (Labusch . 1970f). domain walls in fcrro 



magnets ( Hilzinger and Kronmulierl 1 19761 ) to name just 
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a few. Two bodies in contact form a set of n ran- 
dom contacts per unit area. In the limit of small uni- 
form loads, we can consider a weak-pinning hypothe- 
sis, in which friction results from the fluctuations of the 
forces d ue to individual contacts, and derive a scaling 
theory (fCaroli and Nozieres L 19981: Persson and Tosatti 



119991: ISokolofA |2002at IVolmer and Nattermannl . Il997ft . 
Neglecting for simplicity tensor indices, the displacement 
of a contact at x due to the elastic interactions with the 
other contacts can be estimated as 



J dVG(x-xV(x'), 



(7) 



where the elastic Green function scales as G(x) ~ 
l/(-E|x|), E is the Young modulus (which is typically 
10 — 100 times larger than the yield strength cry), and 
the contact stresses a are randomly distributed with zero 
mean correlations 



(cr(x)cr(x')) = na Y a S(x — x') , 



(8) 



where a is the contact diameter, depending on the 
normal load, the yield stress ay is taken as a mea- 
sure of the stress due to the contacts, and n is the 
contact density. In analogy with other collective pin- 
ning theories, we can estimate the typical length- 
scale L r [usually referre d to a s the "Larkin length" 



( Larkin and Ovchinnikoy . 1979h or "elastic coher ence 
length" (jMiiserl . 12004 iPersson and Tosatti Il999h ] at 
which clastic interaction become important by the condi- 
tion that the typical displacement equals the size of the 
A straightforward calculation yields 



contact (u 



L c = a cxp C 



Ed 
ay a 



(9) 



where C is a numerical constant, and d = ^j\fn is the 
average distance between contacts. The Larkin length 
separates the lengthscalcs for which elastic interaction 
dominates (L < L c ) and the contact interface does not 
deform, from those at which disorder dominates (L > L c ) 
and interface adapts to the pinning-center landscape. For 
two macroscopic bodies in contact, Eq. ^ predicts that 
L c is very large because d^> a and E 3> <Jy , and pinning- 
induced deformations should be absent. One should bear 
in mind that boundary-induced stress gradient can still 
lead to observable elastic deformations, as discussed in 
the next section. 

An important aspect that is missing from the analy- 
sis above is the presence of long-range correlations: most 
contact surfaces are self-affine over several lengthscales. 
This implies that the assumption of uncorrelatcd surface 
stresses made in Eq. ([SJ is not valid. One can, however, 
repeat the Larkin argument for long-range correlated pin- 
ning stresses 



( CT (x)a(x')> = 4 



(10) 



^ v ^ ^ ^ ^ 
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Figure 12 A sketch of the Burridge-Knopoff model. A set 
of factional blocks connected by springs of stiffness Ko are 
attached to a slider moving at velocity v by a set of springs 
of stiffness Kd- 



where 7 is a scaling exponent. Computing the typical 
displacement in this case, we obtain 



E 

L c ■ ot a I 

ay 



(11) 



for 7 < 2, while for 7 > 2 long- range correlations are 
irrelevant ( Fedorenko et all 120061 ) and we recover the 
uncorrelated case in Eq. (j9j). The interesting feature of 
Eq. (fTTj) is that the dependence is not exponential and the 
contact interface can deform even at small length scales. 
The effect of clastic interactions due to the contact of 
self-affine surfaces has been studied analytically and nu- 
merically, revealing that t he solid indeed deform s elas- 
tically due to the contact ( Campana et all 120081 1201 it 
iPerssonl [20061) . 

The role of elasticity in friction is vividly illus- 
trated by the Burridge-Knopoff model for earth quakes 
(|Burridge and Knopofj . ll967l:ICarlson et a/.l . ll994l ) where 
a set of frictional blocks of mass M coupled by springs 
are driven over a substrate (see Fig. |T2j ■ In one dimen- 
sion, the equation of motion for the displacement Ui of 
block i is given by 



MUi = Ko(ui+i + Ui-i - 2ui) + Kd(ui —ia — vt) + /(ii t ) 

(12) 

where Ko and Kd are the stiffnesses of the springs con- 
necting the blocks between themselves and with the load- 
ing plate that moves at constant velocity v. Here f(v) 
is a phenomenological friction force that decreases with 
the velocity of the block and a is the rest length of the 
springs connecting the blocks. 

The Burridge-Knopoff is a deterministic model but dis- 
plays a very rich dynamical behavior with widely dis- 
tributed slip events. The key to the instability is the 
velocity- weakening constitutive law /(w) employed to de- 
scribe the frictional properties of each block. The model 
thus operates at a macroscopic scale and one needs to 
justify microscopically the origin of its constitutive law. 
In this context the model serves to illustrate the concept 
of localized instabilities during friction: on the tectonic 
scale slip is localized on some portions of the fault that 
does not necessarily move coherently as a rigid block. 



C. Mesoscale friction: detachment fronts 

Significant progress in understanding the relationship 
between the dynamics of individual contacts and macro- 
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scopic frictional motion has been achieved with the devel- 
opment of a new real-time visualization method of the net 
area of contact along the entire interfa ce, as discussed in 
a ser i es of papers by Fineberg's group (IBen-David et al 



2010; Ben-David and Fineber 



s group l\a 



ubinstein et al 



20041 . 120061 l2007h . In their experimental apparatus, two 



poly(mcthyl mcthacrylate) (PMMA) blocks were pressed 
together and sheared by a constant force. A similar vi- 
sualization technique has been used for tribological stud- 
i es of a different tra nsparent material, Columbia Resin 
( Nielsen et all I2010T ) . Owing to the transparency of the 
media, it was possible to record the contact area as the 
block were slipping. This method has enabled a number 
of key conclusions on the mechanism of transition from 
static to kinetic friction in macroscopic systems to be 
drawn: (i) the onset of sliding is preceded by a discrete 
sequence of crack-like precursors (collective modes of the 
entire ensemble of asperities); (ii) the transition is gov- 
erned by the interplay between three types of fronts: sub- 
Rayleigh, intersonic and slow fronts; and (iii) a sequence 
of "precursor" events gives rise to a highly inhomoge- 
neous spatial distribution of contacts before the overall 
sliding occurs. The collective behavior of the asperity 
ensemble that composes a frictional interface therefore 
determines the transition mechanism from static to ki- 
netic friction. Imaging the contacts during friction, also 
allowed to record the l ocal values of shear and normal 
stresses. By doing this, IBen-David et al. showed 
that the friction coefficient is not a constant material 
property but it also depends on the way the s ystem is 
loaded locally (jBen-David and Fineberet l201ll) . These 
results call into question many assumptions that have 
ruled friction for centuries: if the onset of sliding occurs 
by the propagation of fronts then elastic deformations 
in the contact interface become relevant and should be 
taken into account by a theory of friction. This may 
even call into question the Amontons law, stating that 
the friction coefficient is independent on the normal load 
and on the apparent contact area. Indeed deviations from 
Amontons' laws have been reported for this experiment. 

The experimental phenom enology is well capt ured by a 
minimal spring-block model (jBraun et "oil. 120091) that de- 
scribes friction at the slider-substrate interfaces in terms 
of rupture and reattachment of surface junctions, which 
represent asperities between rough surfaces. Contrary 
to the above discussed Burridge-Knopoff model of earth- 
quakes where phenome nological friction la ws have been 
introduced, the model of lBraun et al. I (|2009t) explicitly in- 
cludes two most relevant material properties: interfacial 
elasticity and thresholds for yielding/rupture of surface 
junctions. The interfacial elasticity, which was ignored in 
previous commonly used models of friction, defines a new 
velocity scale that is independent of the Rayleigh speed 
and corresponds to slow crack-like fronts ( Nielsen et ali 
\20ld iRubinstein et all 12004 120061) mediating the tran- 
sition from static to kinetic friction. 

The spring-block model (IBraun et al 1 I2009D motivated 
a continuum description of friction between spatially ex- 



tended materials that includes a coupling between the 
bulk elastic deformations and the dissipative dynamics 
at frictional interfaces, which mimic the ruptu r e and 
reattachment of microc ontacts ( Bar Sinai et ali l2012t 
iBouchbinder et all 1201 lh . This promising approach may 
bridge a gap between microscopic and macroscopic scales 
and enable simulations of macroscopic friction processes 
at time and length scales relevant to tribological mea- 
surements. 

The ID sp ring-b lo ck models discussed above 
(IBraun et all 120091: iBurridge and Knopofj . Il967t 
ICarlson et al\ , Il994h have limitations, which do not 
allow a quantitative description of macroscopic friction 
experiments. In particular, the ID models predict an 
exponential decay of elastic interactions within the 
sample, while a 3D description leads to a power law 
decrease of the stress. Nevertheless, recent 2D calcu- 
lation s employing spring-block m odel (|Tr0mborg et all 
l201ll ) and finite-element method ( Kammer et all 20121 ) 
demonstrated that the ID models provide an important 
insight into the mechanism of dry friction between 
spatially extended materials, and allow to investigate 
the effect of system parameters on frictional response. 
It is important to note that, in order to reproduce 
experimentally observed spectrum of detachment fronts, 
the models should incorporate interfacial stiffness and 
local friction laws including b oth velocity we a kening 
and strengthening branches dBar Sinai et 

Mi l2012t 

IBouchbinder et al. [ |20 lit IBraun et all\200$l The mod- 
els assuming local Coul omb- Amontons f r iction at the 
block-substrate interface ( Tr0mborg et aZ.I . l201lh do not 
exhibit slow ru pture fronts li ke those o bserved in various 
experiment s ( Nielsen et all |2010| ; IRubinstein et al] , 
I2004l2006l) . A quantitative comparison to experimental 
data requires 2D calculations with a proper choice of 
local friction laws that is a challenge for future studies. 

Several questions still remain open: how general are 
these experimental results? Do they depend on the spe- 
cific material, or the setup geometry? If not, should 
we revise our general understanding of friction based on 
Amontons' laws? Yet Amontons' laws have worked quite 
well for centuries so they should still be valid at least in an 
average sense or in some conditions. All these questions 
require further theoretical insight and more experiments 
in mesoscalc friction. 



V. A FEW SPECIAL FRICTIONAL PHENOMENA 



A. Electronic friction 



As many theor ists dGoncalves and Kiwi 1999 ; 



Novotnv and Velickv 



1999k 



Persson, 



Persson and Volokitin , 19961: iPlihal and Langreth , 1999 : 



2000b; 



Schaich and Harrisill98lHSokolofj . |2002bl:ISokoloff et al , 
20001 : Volokitin et all 2007 ) have discussed, sliding fric- 
tion over a metallic substrate should elicit electronic 
excitations, giving rise to additional frictional dissipation 
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Figure 13 Temperature variation of the friction coefficient 
across the critical point T c = 9.2 K of Nb. Red squares / 
black dots: tip-sample separation 0.5 nm / several /jm (free 
cantilever). Green solid line: fit by the analytic dependency 
expected from the Bardeen-Cooper-Schrieffer theory. The 
fricti on coefficient is shi fted vertically by 2.5 x 10 -12 kg s _1 . 
From iKisiel et al] (120111 ) . 



besides that due to phonons. Thus for example friction 
on a metal surface should drop when the metal is cooled 
below the superconducting critical temperature, as 
normal electron-hole gapless excitations disappear. A 
first QCM report of this frictional dro p for the sliding o f 
molecular N2 islands on a Pb surface (jDavo et all Il998l ) 
broke the ic e, but also tri ggered a considerable debate 
(jFois et all l2007t iKriml Il999t iRenner et all , Il999t) . 
More recently, electronic friction was demonstrated in 
doped semiconductors, where local carrier concentration 
was controlled through application of forward or reverse 
bias voltages between the AFM tip and the sample in the 
p and the n regions, thus demonstrating the capability 
to electroni c ally t une friction in semiconductor devices 
dPark et all 120061) . 

Very recently, using a pendulum-type AFM probe, a 
clear noncontact friction drop over the surface of Nb 
has been characterized at the superconducting transi- 
tion ( Kisiel et all 1 2 llh . The features observed at this 
transition follow quite closely the predictions by Pers- 
son and collaborators, see Fig. [T3] This ultra-sensitive 
pendulum probe is now ready to be put to work to 
detect the change of electronic friction at other su- 
perconducting transitions. The main progress to be 
expected is now experimental, more than theoretical. 
Promising cases s hould include t he high-Te cup rates 
(iDamascelli et all 120031). or ganics dKanodal 120061) . ful- 
lerides ( Capone et al , 20091: Gunnarsson . |2004) . heavv- 
fermion compoun ds ( Ernst et all l201fjf ) and pnictides 
( Yin et al\ , l2009h . One interesting physical issue that 



these investigations could address is the nature and role 
of strong electron correlations, generally believed to be 
important for superconductivity in most of these materi- 
als. Strongly-correlated superconductivity models gener- 
ally imply that, when not superconducting, the metal 
is not always Fermi-liquid like, therefore with vanish- 
ingly small quasiparticlc strengths and Drude weights. 
Moreover, it should be noted that a recent experimen- 
tal and theoretical study demonstrated the general ex- 
istence at a strongly correlated metal surface of an in- 
sulating "dead layer" with a thickness wh ich for exam- 
ple in V2O3 reaches several nanometers ( Borghi et all 
120091 : iRodolakis et all 120091 ). Another interesting case 
to be studied could be two dimensional metal-insulator 
transitions at some surfaces, such a s the Mott transition 
reported for Sn/Si( l ll) y/ 3 x ^3 (|Modesti et all l2007t 
iProfeta and Tosattil . [20071) . 



B. Magnetic dissipation 

The relationship between nanofriction and magnetism 
at the atomic level is an intriguing side direction. In a 
recent Mag netic Exchange Force Microscopy (MExFM) 
experiment (|Wiesendangerll2009l) . atomic force sensitiv- 
ity to the magnetic state of a surface atom was demon- 
strated for an atomically sharp Fe mag netic tip over 
the ( 001) surface of antiferromagnetic NiO ( Kaiser et all 
I2007D . Besides a different force for the two oppositely 
polarized surface Ni atoms - well explained by the Fe-Ni 
exchange available from elec tronic structure calculations 
( Momida and Oguchil . l2005l) - the results show a surpris- 
ing difference of mechanical dissipation, with a large ex- 
cess of order 15-20 meV per cycle in the antiparallel Fe-Ni 
spin configuration, as compared to the parallel one. The 
Fe-Ni exchange energy is higher in the antiparallel case, 
and the difference can clearly be dissipated by flipping 
the Ni spin. However, direct excitation of surface anti- 
ferromagnetic magnons in the antiparallel tip-Ni case - 
the first obvious possible explanation - is ruled out since, 
owing to strong dipolar anisotropy, the antiferromagnetic 
spin- wave d ispersion of NiO has a gap A ~ 1.5 meV 
- 0.36 THz (jHutchings and Samuelselll972l ) in bulk, and 
one at least as large at the surface. The low- frequency 
oscillatory perturbation exerted by the tip (~ 160 KHz) 
on the surface spin, far smaller than this gap, is com- 
pletely adiabatic, and direct dissipation in the spin-wave 
channel is impossible. Other strictly magnetic dissipa- 
tion mechanisms involving me soscopic s c ale ph enomena, 
such as domain wall motion ( Liu et al\ , Il997l ) , also ap- 
pear inefficient in the atomic scale tip-sample magnetic 
dissipation. One is left with magnetic coupling to surface 
atomic displacements. Acoustic phonons are not gapped, 
both in bulk and at the surface, so they could indeed dis- 
sipate. However, the probe frequency is extremely low; 



and since dissipation vanish es in linear response theory as 
a high power of frequency ( Persson et all Il999l) . a mag- 



netic dissipation mechanism via phonons must involve 
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some hystcretic ph enomena far from linea r response. Re- 
cent theory work ( Pellegrini et all |2010|) suggests that 
the nonlinear response may be related to the attain- 
ment of a strong coupling overda mped spin-phonon state 
very well known in other contexts ( Caldeira and Leggettl . 
Il98ll ). giving rise to an unusual kind of single-spin hys- 
teresis. The tip-surface exchange interaction couples to- 
gether spin and atom coordinates, leading to a spin- 
phonon problem with Caldcira-Lcggctt type dissipation. 
In the overdamped regime, that coupled problem can lead 
to a unique single-spin hysteretic behavior with a large 
spin-dependent dissipation, even down to the very low 
experimental tip oscillation frequencies, just as is seen 
experimentally. A quantum phase transition to an un- 
derdamped regime with a loss of hysteresis and a dra- 
matic drop of magnetic tip dissipation should in prin- 
ciple be found by increasing and tuning the tip-surface 
distance. This experimental check would also help distin- 
guish this interesting spin-phonon mechanism from more 
trivial possibilities, such as a additional dissipation sim- 
ply due to a closer distance Fe tip - NiO surface approach 
in the antiparallel spin configuration. 



C. Carbon nanotube friction 




R(nm) 



Figure 14 (color online). Simulated flow of water (a) inside or 
(b) outside armchair CNTs, and the ensuing friction depen- 
dency on (c) slip velocity and (d) confinement radius R. Fk is 
the kinetic friction force, and A is the friction coefficient. Pan- 
els (c) and (d) include data for graphene s labs, for which 2R 
is the wall-to-wall distance. Adapted from lFalk et al\ (|2010f ). 



Carbon and carbo n nanotubes (CNT ) are widely em- 
ployed in tribology (jChen et all |2003( ). N anotube ap- 
plicat ions are a lso numerous in nanofr iction. iFalvo et al\ 
(|l999h and also lBuldum and Lul (|l999f ) discussed the pos- 
sibility to slide or roll nanotubes on a surface. Recently, a 
large longitudinal-transverse anisotropy of AFM friction 
on surface deposited CNTs has been explained precisely 
by the contrast between the longitudinal tube rigidity, 
against t he transverse softn ess manifested by "hindered 
rolling" (|Lucas et all l2009f ) . 

Zettl and co llaborators ( Cumings and Zettl . l2000t 
iKis et all I2006T) demonstrated ultralow friction expe- 
rienced by coaxially sliding multiwall CNTs. Coax- 
i al CNT sliding also i nspired num erous simulations 
(jServantie and Gaspardl 120031 l2006alfbh . More simu- 
lation identified curious - even if at present rather 
academic - frictional peaks at selected sliding speeds 
( Tangnev et all 120061) corresponding to the parametric 
excitation of quanti zed nanotube breat hing modes. Re- 
cent theory work ( Zhang et all l2009h also discovered 
that the twofold degeneracy of the breathing modes can 
cause chiral symmetry to break dynamically at these fric- 
tional peaks, so that even purely longitudinal coaxial slid- 
ing of nonchiral tubes can generate angular momentum. 
Other mechanical /rheologi cal propertie s of C NTs have 
also been probed by AFM (jPalaci et all |2005[) . 

Water wetting and the interfacial friction of water 
in CNTs have also been studied for various purposes. 
The observation of flow-generated voltages growing loga- 
rithm ically with velocity of ion-rich water ( Ghosh et all 
I2003D appears to be a manifestation of electronic friction 
which has found various explanations including stick-slip 



of ions embedded in high viscosity water near the tube 
(Pcr sson et al ., 2004) or a statistical consequence of the 
flow-induced asymmetry in the correlation of the ions, in 
the a mbient fluid as seen by the charge carriers in the 
CNT (jGhosh et aZ.1 . 12004^1 . 

Recent studies have focused on disentangling confine- 
ment and curvature effects on water friction inside and 
outside CNTs, showing that the friction coefficient ex- 
hibits a strong curvature dependence. While for a flat 
graphene slab friction is independent of confinement, it 
decreases with CNT radius for water inside, but increases 
for water outside, see Fig. [141 The friction coefficient is 
found to vanish below a threshold diameter for armchair 
CNTs. A "superlubric" , structural origin of this cur- 
vature dependence associated with a curvature-induced 
incommensurability betwe en the water and carbon struc- 
tures, has been proposed ( Falk et aZl . l2010h . 



D. Friction in colloidal systems 

Handling matter with static periodic fields generated 
by interfering lasers adds a new exciting actor to the 
realm of toy systems that display real physics. Trap- 
ping and handling colloidal particles with intense photon 
fields offers the possibility to change parameters freely, to 
compare directly experiment with theory, to test theoret- 
ical predictions, and to visualize directly in simple cases 
the intimate mechanics of sliding friction ( Vanossi et all 
120121 ; IVanossi and Tosattl 120121 ). In the novel approach 
freshly inaugurated by Bechinger's group ( Bohlein et all 
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|2012| ). a 2D close- packed crystal of charged colloidal par- 
ticles is forced to slide by Stokes forces in the presence 
of a laser-generated static potential. Unlike conventional 
sliding crystal surfaces, two-dimensional lattices with dif- 
ferent symmetries, lattice spacings and corrugation am- 
plitudes can be constructed at will, realizing for example 
commensurate or incommensurate matchings, quasicrys- 
tal substrate geometries, and possibly "disordered" ge- 
ometries too. 

Besides and above all that, colloidal friction provides 
an unprcccdent real-time visualization of the full fric- 
tional dynamics at play. Unlike AFM, SFA ad QCM 
which provide crucial, but averaged frictional data such 
as the overall static and kinetic friction, mean velocities, 
slip times, etc., the colloidal experiments photograph the 
actual time-dependent motion of every individual parti- 
cle during sliding, an exquisite privilege restricted so far 
to the ideal world of MD simulations. Transitions be- 
tween different dynamical states become experimentally 
accessible and can be analyzed and related to the de- 
tailed particle motion. These and other goodies lie in 
the future. 



VI. CONCLUSIONS 

Among provisional conclusions of this Colloquium we 
mention: 

(i) Despite limitations, all levels of modeling and sim- 
ulation are highly informative and predictive. 

(ii) The PT model is, despite its deceiving simplic- 
ity, extremely useful in understanding many aspects of 
nanofriction. Its extreme success tends however to hide 
the actual complexity of the phenomenon. 

(in) Multi-contact models are instrumental in describ- 
ing mcsoscopic friction and fracture. One main problem 
of these models is the multiplicity of empirical parame- 
ters they involve. Other open problem are under what 
conditions these models yield realistic stick-slip friction 
as opposed to smooth sliding, and the exact role of the 
elastic interactions between contacts. 

(iv) Molecular-dynamics simulations are good and in- 
formative for qualitative descriptions of atomic stick-slip. 
Open problems are the high speed of observed change 
from stick-slip to smooth sliding; the potential artifacts 
introduced by unrealistic dissipation mechanisms of Joule 
heat, and more importantly the simulation size and time 
limitations, in particular the complete omission of slow, 
logarithmic relaxations and ageing. 

(v) Prospective mechanisms for the control of friction, 
such as mechanical oscillations or a phase transition in 
the substrate, are suggested by model studies, and are 
presently under experimental scrutiny. 

It is worth mentioning in closing that there remain 
fully open problems at the very basic theoretical level. 
First of all, we do not have a proper theory of friction, 
namely a theory where the frictional work could be cal- 
culated quantitatively (not just simulated) in all cases - 



they are the majority - where linear-response theory is 
inapplicable. Second, while MD simulations can to some 
extent be used in lieu of theory, they have been so far 
strictly classical. Future work should include quantum 
effects and gauge their importance. 

There are many more outstanding challenges left in 
nanofriction. Among them: 

• Bridging the gap between nano and meso (macro) 
scale friction: multicontact systems. 

• Mechanical control of friction in multicontact sys- 
tems. 

• The ageing of surface contacts at the nano and 
macroscales. 

• Role of wear and adhesion at the nanoscale. 

• Rolling nanofriction: besides the known case of 
nanotubes: does it exist, and how to distinguish 
between rolling and sliding? 

• Friction in dislocations and in granular systems. 

• Friction in biological systems (motor proteins, cells 
membranes and pores, etc.). 

Lively progress along these and newer lines is to be ex- 
pected in the near future. 
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